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A3 STBACT 

In this thesis certain aspects of radio wave 
propagation through the ionosphere are investigated with 
the help of ray tracing using a digits,! computer for the 
purpose. A computerised three dimensional ray tracing 
program is used Toased upon the Haselgrove equations and 
employing the Appleton-Hartree equation for the index 
of refraction, j The program is based on the work of 
Jones (1966). This routine has been used to predict the 
effect of elevation angle and frequency of transmission 

upon ground distance, maximum height of reflection, 

% 

various path lengths and ray path trajectories of HI’ 
and VHP radio waves through regular and disturbed 
ionosphere with and without considering the effect of 
imposed magnetic field. Effects of collision are not 
presented in this paper. A North-South propagation 
path along constant geomagnetic longitude is considered. 
The transmitter is located at Kanpur (26.5°N, 80,5°E) 
and receiver is either on ground or on a satellite at 
height of 1000 Kin above ground.^ 

Suitable models for earth’s magnetic field, 
ionospheric electron density and perturbation to the 

ionospheric electron density have been used.'' '”^ | 

' ' ' . ' ' 1 

Gaussian perturbation (torus) to the ionospheric electron | 

' ■ ' '' X I 

density model has been introduced for this work. '‘The 



basic fea.tures of the ionospheric irregularities 
and phase path characteristics for oblique ray paths 
are reviewed along with the xmderlying theory of the 
ionospheric ray tracing algor it hms . v 

The result of this simulation indicate that in 
general: 

(i) Ray path trajectories are significantly affected 
in the vicinity of torus, the rayswhose wave 
reversal occur near the center of the torus are 
the most significantly affected. 

(ii) Ray path symmetry about the point of reflection’ 
is not disturbed appreciably. 

(iii) There is significant change in various path 

lengths and skip distance. : 

(iv) The highly non-concentric isoelectron contours ! 
caused by the presence of the torus may result 

in significant large scale focussing/defocussing 

i 

of the radio waves. In general, ray reflected j 

! 

from the center portion of the torus are i 

focussed. j 

(v) For any undisturbed ionospheric electron density 
profile and signal transmission frequency above * 

' i\ 

critical frequency, a cusp is observed in a | 
curve depicting voxiation of phase path length > 
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and ground range with elevation an^e , whereas 
a nose appears in a curve depicting variation 
of group path length and ground range with 
elevation angle. In former case the slope of 
the curve changes gradually whereas in the 
latter case the slope of the curve changes 
abruptly . i 



Chapter 1 


lETTRODUCTION 


1.1 General ^ 

Ray tracing is a technique of finding the path 
followed by radio waves from transmitter to receiver 
through an ionized medium in, terras of geometrical 
optics. For ionospheric problems, ray tracing can be 
used to determine such basic parameters as the 
direction of energy flow, ground range, maximum 
height of reflection, phase path, group path, actual 
path length etc. From those basic data, many other 
desired quantities can then be deduced. A three 
dimensional ray tracing program developed by R.M. Jones 
of ESSA (1968) written in CDC 3600 PORTRAY, which is 
similar to FORTRAN IV is used to trace radio rays 
throu^ the ionosphere. This could be converted to 
other applications by substituting an appropriate 
index of refraction subroutine. 

Since the first years of ionospheric stiidios, 
as early as 1920's, certain forms of ray tracing have 
been used for ionospheric problems which involved 
laborious calculations so as to get an idea of 
refraction according to the laws of geometrical optics 
but within tlae past decade the increasing availability 
of high-speed digital computers has greatly stimulated 
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the development and application of ray tracing process 

"by permitting extensive nse of voluminous and complex 
calculations which earlier could he done only in rare 
cases and at %he expense of much time and labour. 

iTesd for Ray Tracing 

To determine the amplitude, polarization, 
relative pha.se, time of flight, absorption, Doppler 
frequency shift, skip distance, maximum usable fre^ •. 
quency (MUF), etc. of a radio wave, it is necessary 
to determine the ray path between transmitter and 
receiver. At first glance, it might appear that it is 
easy to determine the path along which the energy 
flows. However, this is not tho case in practice. 

Hor a few simple ionospheric electron density lT(h) 
profiles, in which electron collisions, the earth's 
magnetic field .and ionospheric irregularities are 
neglected, it is possible to determine ray path 
analytically. In general, however, the paths have 
to be determined by a step-by-step process. There 
are varying degrees of sophistication in ray tracing. 
In simplest form of ray tracing the ionosphere is 
considered as a plane mirror, and ray paths are made 
up of linear segments. In more sophisticated ray 
tracing programs the effects of electron-density 
variations (sufficiently slow for ray theory to 
hold), magnetic-field changes, collision frequency, 
perturbation to the electron deisity are all taken 
into account. 



3 


1.3 Scope 

Ray tracing techniques can he accurately and 
efficiently Cipployed for studying radio ray path 
trajectories in ionosphere and troposphere at HP hand, 
ionospheric phase distortion at UHP hand and other 
effects such as absorption, Doppler frequency shift, 
hearing deTieition, error in angle of transmission, 
etc. This work makes use of ray tracing techniques 
to study such parameters a,8 ground range, phase path, 
group path, etc. for various elevation angles and 
frequencies of transmission in HP and VHP hand. 

Stu^’Pes have heen made hoth with and without the 
effect of imposed magnetic field of the earth. The 
effect of collision has not heen considered. Howevor, 
the same can he included hy using suitable collision 
frequency model to estimate absorption. An analytical 
model of electron density has heen used to represent 
the electron density contours. A perturbation model 
is used to represent the effect of torus upon the 
electron density contours. Since the primary 
objective of introducing perturbation is to deduce 
order-of-magnituri e characteristics of the ray path 
trajectories, the presence of the earth's magnetic 
field has been neglected while studying the effect 
of perturbation. Ray path trajectories and phase 
path characteristics for oblique angle rays have been 
plotted by IBM 1800 computer to depict effect of 
elevation angle and frequency of transmission through 
a regular as well as perturbed ionosphere. 
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, Chapter 2 

IONOSPHERE 

■5 

2»1 Ionospheric Structure 

The ionosphere is that region of the earth's 
atmosphere l 3 ring approximately between 50 Em and one 
earth radins (6370 Ere), in which sufficient ioniza- 
tion exists to enfluence the propagation of radio 
waves. The ionosphere is sub-divided into D,E and 
P regions within which layers of electrons exist .The 
chief characteristic of each layer is the critical 
frequency which is the highest freniaency reflected 
by the layer for vertical incidence. It depends 
upon the electron densitjr. 

2.1*1 Regions and I'lyers; 

The approximate locations of regions and the 
layers that may exist within them are given in Table 
2.1 below. 


TABLE 2.1 


Region Height 

Layer 

Approxims-te 

Approximate Day 


Era 


height, Em 

time Electron 
Density i2i-3 



C 

€5 

10® 

D .. 

50-90 

L 

75-80 

0 ^ 1 
0 

H 1 

E 

90 to 

El 

110 

H 

H 

0 . 

H 


120-140' 


100' 

Unpr e'd i ct 'abl e ' 


120-140 

Fl,Fli 

200 

2x10^^ 

F 

and 


10^^ 


above 

F2 

above 250 
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C-laycr (50-70 Pin) is produced by cosnic rays, 
vibereas the other layer, are prodiiced by solar radia- 
tions, with the exception of sporadic S layer (E„) 

D 

which is produced by meteoric particles. D-layer 
exists in daytime with its maximum electron density 
nea.r 80 Zm with an average '.thickness of 10 Km. In 
D-layer degree of ionization depends on the altitude 
of sun and therefore it disappears at night. El-layer 
is always present by day in all seasons and in all 
geographical locations and its behaviour is subject 
to close solar control. The sporadic E-layer (S^) 
consist largely o'-' relatively dense patches of 
electrons of the order of a few tens of kilometers 
in horizontal extent. The critical frequency of 
E-laj^-er is given by 

f^E = 0.9[(180 + 1.44R) cos (2.1) 

where R is Zurich sunspot number (SSE) and X is 
Zenith angle. 

Fl-layer (200-250 Ktd) with day time average 
thickness of 20 Km', combines with E2-layor at night. 

It undergoes seasonal as well as solar cycle varia- 
tions. It is more pronounced during summer than 
during winter and at high sunspot numiber. The 
critical frequency of El- layer is given by 

f^El = (4.5 + O.OIR) cos^'^X (2.2) 




To a first approximation the exponentv. in equation (2.2) 
varies from place to place and with season. It is less 
during sunspot minimum than at sunspot maximum. P2- 
layer (200-400 Em) is hy far the most important 
reflecting medium for HE radio waves. At night it 
falls to a height of 300 Em and comhines with Fl- 
layer. Its height and ionization density vary largely, 
and depend on sunspot cycle and ijinlergo diurnal and 
seasonal variations. Its critical frequency varies 
linearly with SSN (for R = 0 to 150) .and thereafter 
it flattens out. It is given by 

f^F2 = 5.4(1 + 0.0064R) (2.3) 

2.1.2 A Chapman layer; 

E and F layers:. are not very distinct and with 
sufficient accuracy it can be assumed that during 
day li^t hours E-layer is a Chapman layer and obeys 
the law as given below. 

N = Nq exp ■ 2 -[l-z-sec xe3cp(-z)] (2,4) 

where 



X = Zenith angle 
H = scale height 

h = height where electron density is computed in 
relation (2,4) 
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= height of m axisaijm ion proinction when 
SLin is overhead. Reference height. 

The critical frequency of E-las^-er (with H = 10 Km, 

Rq = 115 Km, Nq = 2.8x10^ (m“^) is 4. 7( cos X MLz 

during day and is of the order of 0.5 MHz at night. 

2.2 Ray Paths in the lonosphc-rc- 

A radio wave tre,velling from transmitter (T) 
to receiver (R) through ionosphere oJ-ong a curved 
path is depicted in Figure 2.1. The basic expres- 
sions for phase path (P), group path (I”) and ground 
range ( D ) are 

R 

P = ^ y, ds 

R 

P* = / p’ ds (2,5) 

T 

D H„a6 

where p = phase refractive index of the ionosphere 
p'= group refractive index of the ionosphere 

= “l^(pf) 

- u + f 

- p + I df 

d = an element of path length and 
s 

d e = angle subtended at the centre of the , ,, 
earth by ds as shown in Figure 2.2. 
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Bach integral is ^taken over the entire path between 
transmitter and receiver. In absence of magnetic 
field and electron collision, the group refractive 
index is inverse of phase refractive index. 

5 

Applying Bouger’s rule to the geometry of 
Figure 2.1, we have 


\XT sin i = Rq cos (2.6) 

From the ray geometry of Figure 2.2, the elements of 
angle -and length are related by the equations 


ric! - ^ ^ Q _ 

sin i cos i 


(2.7) 


and therefore, from equations (2.5), (2.6) and (2.7) 
we get 


P = [i ds 
T 

= f _ii_dr_ 

j cos i 

^V 2 , 

= 2 I I". mi; r 

E (n= cos= 

0 ^ 

P* = /^ ds 

" T ' ' ' ' ■ 

= f M 

^v 

1 r dr 

- Rq cos^ 


(2.8a) 


i 


(2.8b) 



2 
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R 

D = ^ Ho dO 


where is the value of r at the vertex of the 
path . 

2.3 Ionospheric Irreigularities 16>-18 

The ionosphere has been under close observa- 
tion for over 35 years. Until 1953, investigations 
have been concerned mostly with the study of long- | 

period variations such as diurnal and seasonal i 

changes in certain parameters, in particular, the | 

height and maximum ionization of the main regions | 

(E,F1,P2). During past feTir years, with the advent I 

of artificial satellites, it has become possible to j 

I 

scan a largo section of the ionosphere in a few | 

minutes. This greatly simplifies the study of | 

ionospheric irregularities and much useful infonria- i 

tion has already been obtained from observations i 

! 

of satellite signals. j 

] 

[ 

2.3.1 Basic Features t 

Experimental observations indicate tha,t the I 

ionospheric irregularities are a permanent features : 

of the day li^t-hoprs ionosphere. The ionization I 


2 / 

R_ 


Rq cos A dr 




- r: 


CO£ 


A)- 


(2.8c) 
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onhancemont is to extend thi- ough. the entire 

F2 region and even into D and E layers. The irre- 
gularities have an electron density which is 
typically 10 percent higher than that of its 
h a ckground . 

2.3.2 General Characteristics; 

Evidence for the existence of much larger 
irregularities in the E region, with dimensions of 
several hundred kilometers, was given hy Mimro 
( 1950 , 1958 ) by comparing ionospheric virtual height 
records taken at different times and different 
places. Similar measurements have been mad. e 13 ^ 
several workers by observing changes in the time of 
travel, the amplitude, or the frequency of radio 
signals received from two or more sta,tions . These 
observations only reveal the presence of irregulari- 
ties that are m.oving horizontally with a mean speed 
of 500 Km/min. The irregulariMes generally have 
horizontal dimensions of a few hundred kilometers. 

A few isolated irregularities as large as 2000 Em 
are observed by Chan and Villard (1962). The 
presence of large scale irregularities in the 
E region has also been shown by Vitkevitch and 
Eokurin (1958) and by Lawrence et al ( 196 I) from 
observations of slow, irregular variations in the 
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apparant positions of radio stars. The horizontal 
dimensions of the irregularity were generally between 
200 and 400 Km with electron density variations of 
about one percent. 

The presence of large irregularities in the 
ionosphere was also dodiiced by Titheridge (1963) 
from calculations of the number of free electrons 
between an observer and a satellite, based on 
measuroments of amplitude of signals received from 
a satellite Explorer 7, recorded at Ardmore 
(lattitude 37.01° S, longitude 174.98° E) 17 miles 
Southeast of Auckland from Nov I960 to Aug 1961. 

He found tho,t there was no seasonal or 
diurnal variations in their occurrence a.part from 
a possible tendency for feweh to be observed in the 
early morning. Irregularities can occur singly or 
in groups of similar size or small ones may be 
superimposed on large ones. The sizes of those 
observed ranged from 5 to 500 Km. 

The total electron content of the irregulari- 
ties^ measured along the line from satellite to the 

-] X 15 

receiver varied from about 8x10 to 15x10 electrons 
per meter square in winter. The seasonal change was 
caused by a decrease in the total electron content 
of the ionosphere from about 2x10^^^ electrons/m ^ in 
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17 / 2 

sujtiEier to 1x10 electrons /in in i^inter. In botla sea- 

sons then, the electron contents of the irregulari- 
ties varied from about 0.04 to 7 percent of the 
total electron content of the ionosphere. 

Titheridge (1968) found, by comparing the 
nimber of irregularities observed using satellites 
a,t different heights, tha,t the number of irregulari- 
ties per unit volume at any height ms approximately 
proportional to the density of the background ioni- 
zation. This suggested that the percentage fluctua- 
tions in the electron density are approximately the 
same at all heights above 200 Km. The irregularities 
were found to exist at heights ranging from 180 to 
750 Km. The horizontal sizes varied from 75 to 
520 Km. The denser irregularities occurred only 
near the peak of the ionosphere between 250-400 Km. 

The irregularities lay on a straight line tilted 
down to the south at an angle of 15 degrees. 

2.3.3 Shape ; 

Titheridge (1963), using the effect of the 
Zenith angle of the satellite on the calculated 
content and size of the irregularities, determined 
their approximate shape. The irregularities were 
not elongated in the direction of the magnetic field 
but occurred in the form of horizontal slabs with a 
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vertical tkickness less than one third of the 
horizontal dimensions, A typical irregularity has 
the shape of an oblate spheroid. 

2.3.4 Direction of Motion®; 

Investigations into the phenomenon of 
travelling ionospheric disturbances (TID) in the 
P-region made by Heisler (1963) at the Radio 
Research Laboratories, University of Sydney, 
Australia, reveals that small perturbations travel 
■'•J'ith slower apparent speeds than large perturbation 
and possess a dominant East-West component of move- 
ment, Their North-South movement is very small. 

By contrast large perturbations have a dominront 
North-South Seasonal component. The difference 
may not be related entirely to perturbation size 
but could be the result of ambiguities and limita- 
tions in observational techniques. It was suggested 
that the estimated speed of a perturbation could be 
a phase speed, group speed or combination of both. 

The large perturbations constitute a major 
ionospheric effect. They are quite prevalent parti- 
cularly during day light-hours, and in winter months 

cause almost continuoiis changes in values of f P2 

0 

(critical frequency of P2 layer). Their movement 
is primarily horizontal, and the associa,tod vertical 
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component of progression which is always observed 
is almost certainly due to a forward tilt in the 
wavefront of the disturbance. 

It seems probable that there 3.re two 
components of motion of perturbations in the 
P-region, an irregular North-South seasonal compo- 
nent which is Diipor imposed on a steady East-West 
drift of ionization which changes diurnally. There 
is considerable evidence to suggest that the 
latter is probably due to North-South electro- 
static fields transferred from the dynamo current 
foci. North-South component of ionization in 
E region is constrainod by the earth's magnetic 
field and can orilj take place along the field 
lines, except by application of an external electrO' 
static field. The observed phase and group veloci- 
ties suggest a wave structure and speed of propaga- 
tion, through medium, of 160 m/s. 

The la,rge ionospheric irregularities 
generally occur in series of throe or moro of 
similar size and shape. This suggests that they 
are caused by a wave of disturbance propagated 
through ionosphere, with a wavelength of about 
100 Km. This was related to the travelling distur- 
bance observed by Munro (1950,1958), which had a 
Spatial length of a few hundred kilometers and were 
attributed to the travelling pressure waves in the 
ionosphere. 
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Cliapter 3 

pha.se path chapacteristics for oblique .ray paths 

3.1 Introduction 

Phase path is a key parameter in theoretical 
discussions of a nujnber of practical problems, 
particiilarly those involving Doppler freauency 
shifts loroduced bj:^ relative motion betxj'een sender 
and receiver or b 3 r changes in the ionosphere arising, 
for example, from Travelling Ionospheric Disturbances. 
For a typical HF oblique ray path from a transmitter 
T to a point of reception R it is possible to calcu- 
late the length, P, of phase path. Fox- a specified 
ionospheric electron density model P can be found by 
numerical integration in a sirLtable computer ray 
tracing program or, for some simple models, by 
evaluation of an appropriate analytical expression- 

3-2 Relation between Phase Path Length and Ground 

5 

Range 

Phase path (P) is a function of frequency (f) 
and elevation angle ( A ) and can be written as 

P = P(f, A ) (5.1) 

It is assumed that the ground range D between trans- 
mitter and receiver remains constant. ifhen the 


frequency is changed by df, the angle of arrival A 


f 
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meas-ured from the horizontal, miist therefore change 
hy an amount* d A to correspond to the now ray path. 

Total differential of relation {3.1) can he 
written as 

dP ' IP [-^ oi 

df = 3f 3A df 


The total change in phase path for a given chanfie 
in frequency is given hy equation (3.2) subject to 
the condition that is zero (where D is the 

groiin'" .range). Since D is also a function of frequ.ency 
and elevation angle, we can write, similarly 


dD 3 D ^ 3 D d A 

3 ? = 3f + 3X 


(3.3) 


this gives 

dA ■ 3 D / 3 D 

^ = “ 3T/ U 

¥e then have 

dP _ 3 P 3 D/ 3 P / 3 D>. /, 

Tf “ ~ 3 f^ 3A ^ 3A ^ 

using equations (2.8a), (2.8b) and (2. 8c), we get 

P-D cosi4 = 2 / ^(u^ - Rq cos^ A )^ dr 

^0 (3.5) 

Partial differentiation of equation (3.5) with respect 
to elevation angle, gives 

^v 

- cos A ^ + D sin A = 2 - — / i(p^r^-RQ cos^ A )2 dr 
3A 3A Rq 


3P 

3A 
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Here both, the integrand and the upper limit of the 
integral are functions of A and therefore the 
differentiation must be C8,rried out according to 
standard rule, i.e., that if 


• 13(7) 

f(y) = / ^( 3 ^, 7 ) dx 

a(y) 


then 


3 y 


a(y) 


^F(x.y) 

3 y 


dx-P(a,y) ^ + P(b,y) 


In our problem, we have ^ = 0 since the lower limit 
is constant, and P{b,y) = 0 because p r.^ = cos A 
at the apex. Hence 


3P 


cos A 



+ D sin A = 





Rq cos a sin A dr 
rCp'r^-R^ cos^a"^ 


= D sin A 


and 


therefore, we 


3_R 
A 


have 



(3.6) 


For a fixed frequency the ground range D and 
phase path P are related through the elevation angle 
as shown in equation (3.6) above. The slope of the 
plot for a fixed frequency f, with D as abscissa and 
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P as ordinate, is cos A . This is true for any 
profile. 

If 'f' is less than or equal, to the vertical 
ci’itical freqiiency f , P ten.ds to twice the phase 
height as elevation mgle approaches 90 deg and 
groun^’’ range approaches zero. The tangent to the 
oi.’rve is horizontal at D = 0. If f>f.,, there is 
a uinium value of ground range (the skip distance) 
at which, from equation (3.6), P is also minimum. 


3 . 3 Relation bet w een P hase Path Len gth and G-r oup 
Pa t h Length ^ 

Prom equations (3.4) and (3.6), wo get 


dP ■ 3P '■ 
I? = rf - 


cos 



(3.7) 


Partial differentiation of equation (3.5) w.r.t, 

'f', holding similar arguments to those in paragraph 
2, gives 

^v 

1^ - cos = 2 / -|p[ ^(p^r^-R^ cos^ A )'^]dr 

"^o 

Assuming that the ionosphere is spherically 
stratified about the centre of the earth, i.e, that 
the electron density is a function of only r, the 
distance from the center of the earth, and neglecting 
earth’s magnetic field and electron collisions, the 
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ref ract i'^T'e index at diotiiiice i* froni the center of 
the earth can be written as 


Hence 


= [1 - 


f^ 




There fore, 

■ .A 


3f 


COS^ 


£D _ 2 r( 1 - ) dr 

9f f ^ (ji^r^-R^ cos^A)^ 


(3.7a) 


xr- - -^O 

Using equations (2.8a), (2.8'b), (3.7), and (3.7a), 
we get 


- 5(^'’ ~ U) (3.8) 

■Above relation holds good for any profile. 

For props.gation via the ionosphere P' is 

always greater than the geometrical path length s 

between transmitter and receiver whereas P is less 

than s. Hence from (3.8) the slope of the curve is 

everywhere positive. There will be cusps in the 

curve wherever^ becomes infinite. Since 41 

df ai 

becomes very large as P’ becomes very large, a 
vertical line through f^^ must be tangential to the 
(P,f) curve; this result applies for any ground range. 
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^ ^ Relation between Group Path LeriA^th and Ground Ran,g:e ‘^ 

For a plane earth, and horizontally stratified 
ionosphere, Breit and Tuve’s theorem gives the relation 
between the effective (group) ohlique-path length P' 
travelled hy a radio-wave packet and the corresponding 
ground range D as 

P ’ = D sec A 

(3.9) 

or D = P ' cJ5s A 

where ^ is the angle of elevation of the ray path 

relative to the ground. Provided that the electron 
density within the layer is a, function of height only 
and that there is no imposed magnetic field, this 
relation is rigorously true for any electron density 
profile. 

Vfith a curved earth the relation between P^ 
and D is no longer independent of profile and the 
simple form of the foregoing equation cannot he used,. 

It is possible to place close limits on any one of the 
quantities P*, P and A when the other two are known, 
■under the same assumption that the earth's magnetic 
field is neglected. We assume the equivalent of 
horizontal stratification, that is, that the electron 
density is a function of r only, where r is the radial 
distance from the earth’s center. 
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For the salce of definiteness, we consider 
the limits that can he put on P‘, given D and A . 
Consider the ra^*" po,tli geometry of Figure 3.1 whore 
the parts of the path below ionosphere, denoted by 
a stibscript 2, are related b 3 r the rigorous equation^ 

p 2 = 2RQSin(D2/2R^) sec( A + D2/2 Fq) 

The limits on P£, the effective (group) path within 
the ionosphere, are obtained as follows. 


Pj_ = /p' ds = I Y 

Also 

ds _ r d 9 
^ M- sin i 


Therefore 




2 T 

r d 0 

+ h^) sin io 


since 

r sin i = (Rq + h^) sin i^ 


All the abve integrals are taken over the path 
"Within the layer Now suppose that we can put 
definite limits on r within the layer, such that 


•mxn 


’max 


It follows that 

r^ . D-, 

min 1 


^c(^o+ ^0 


p » < 

J^l - 


r" D, 
max 1 


R^(Rn + h^) Sin i 


o‘ 


'o 
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since RqCosA = (R^ + h^) sin i^ 

■■re hnve, finally 

2 2 

r . T 

D sec A < P * < — Pn sec A (3.10) 

t 32 1 - 1 — -0 2 1 

j1 s iVq 

ThriS Breit and Tuve’s theorem ma3?' he written in the 
form. 

= a secA (3.11) 

where a = 1 for a plane earth and slightly greater than 
1 for a curved earth. 

If we assamie that the base of ionosphvore extends 
down to the earth then 



For a path involving reflection at the ionsphere, P 
cannot be less than 2hQ,' where h^ is the height of the 
base of the ionosphere, but is less than s. 



25 


For a finite groimd range P is therefore always 

finite, ani mist approach sone finite limiting value P^^ 

for the penetration condition at which P' becomes 

infinite. An infinite value of P' for a finite value 

of D can occujr oni.y at A = 90^, as nay be seen from the 

conditions (3,11) which apply for any profile. From 

e q.uat io ns (2.8a) an-’’ (2.6c), we get 

^v ( |j,^ r^ - cos^ A )'‘^ 

P - D cos A = 2 / dr 

Ro r 

Hence, when A = 90*^ 

P^ = limit [2 / p dr] = 2h^ 

r^-^ (Rq+ la^) ^o 

where h is the limiting valu.e of phase height, i.e., 

O 

the pho.se height for f = f^. The value of P^ is 
therefore independent of the ground range. 
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Chapter 4 

AITALYSIS TECHNIQUE 

The efff-ct of elevation angle and frequency 
on HE radio transmission through a regular as well 
as perturbed ionosphere may he accurately and 
conveniently simulated through the use of digital 
ray tracing techniques in which the effect of electron 
density variations (provided they are sufficiently 
slow for ray theory to hold), magnetic field changes 
and electron collisions are taken into account* A 
three dimensional ray tracing computer program based 
.upon work by Jones (1966) has been used in this thesis 
work. The program utilizes the ray tracing differen- 
tial equa.tions developed by Haselgcove (1954) and the 
Apple ton-Har tree equation for the index of refraction 
in the ionosphere. In addition to determining the 
three dimensional ray path, program was also us^ to 
determine the effects of ionospheric irregularities 
upon the ray path trajectories, ground range, phase 
path, group path etc. Phase path characteristics 
for oblique ray path have also been studied. 

4.1 Ray Tra cing Relationships 

The path followed by a radio wave through an 
inhomogeneous medium, such as the ionosphere, is 
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described by a set of differential oqTjic,tions relating 
the coordinates of a point on the ray path. Haselgrove 

has derived these differential conations (for a 
general coordinate system) from Hamiltonian optics 
using Permats principle of stationary ' time . 

Fermat principle states that a ray travels 
between two points along a path such that the time 
taken is minimized. This principle is stator? in 
the notation of the calculus of vari'rtions .ass 

6 / M ds = 0 (4.1) 

iirhere M = the ray refra.ctivo index 
= p cos a 

p = the phase refractive index 
a = the angle- botvreen v^avc normal and 
wave direction 

ds = an element of path 1.. ngth along the 
ray path. 

Applying Euler's varirtjonal principl.e to equa- 
tion (4-1) expressions for the curvature of the 
ray may be obtained in terms of M and its deriva- 
tives with respect to position and direction. But 
M can be found only when the wave normal direction 
is known; hence the latter must be included as an 
extra set of variables. However, those- equations for 
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the curvature of the ray mp.y he transformed into 
ones uhich give the rate of cho.nge of wave normal 
direction in terms of the-phase r'-.-f roctive index 
p.nd its gradients. The rpy refractive index a,nd 
ray direction are thus eliminated from the ray 
tra. cing equc.t j ons . 

The Haselgrove equations in sphi^rical coordi- 
nates r , 6 > f are ; 


in 

dt 


dt 


d£ 

dt 


p p* 
1 


p p' r 
1 


[V - ii-iH ] 
r 8 v 

r 


-[Vg ] 


p p'Tsine 3V, 


s-Ve 

[V ] 


<p 


(5.2) 


dV^ 

dt" 

dt 


1 ?5u ^0 am 

jL. + y — + Vq 3 sin e ^ 

p» dr B dt ^ 


dt 


Ir 1 iH. 

r'- p' 98 


- Vast + \ 


'dt 


9 


dY(p 

~ 


•[ 1- iii _ sin8 If “ 3^ cose H-] 


r sine 0 


9 


9 


dt- 


where 

p = real ppytk of the phase refractive index 


of the medium 
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VI* = Tt^.l part of the group refra.ctive index 
of the medium 

’ = p + f 

9 f 

V = a vector, directed normal to the phase- 
front of the ray, of magnitude p with 
Vg and its components in the r, 

0 and 9 direction, rL-spectively 

t = group path length (the independent va.riable) 
f = frequency. 

The set of first-order partial differential equations 
are solved numerica,lly to give the locus of the ra,y 
as it traverses the medium which is characterized 
by its parameters such as phase refractive index 
<and group refra,ctivo index. The ray tracing 
procedure is started with initial conditions (i.e,, 
height of transmitter, azimuth angle of transmission, 
and elevation ongle of transmission) and subsequently 
the differential equations are solved numerically 
using ran Adams-Moulton method with Runge-Eutta 
starter. 

In addition to the six basic equations 
necessary to calcula,te the ray pa.th, the program 
adso simultaneously integrates several other 
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differential equntions that are used to determine 
the Doppler frequency shift, and the absorption that 
the radio wave experiences as it propagates through 
the ionosphere. The relationships used to detennine 
these parameters are given in Sections 4.6 and 4.7 
respectively. 

4.1.1 Refractive Index of the Ionosphere ; 

Prediction of ionospheric effects upon 
radio waves requires the knowledge of the complex 
phase and group refractive index of the medium. 

Since 

r,* = p + f (4.3) 

both op and p ’ are known if p is available as a 
function of frequency. 

In the presence of the earth's magnetic 
field and electron collisions, the complex phase 
refractive index of ISie ionosphere is given by the 
Appleton-Hartree equation 

p ^ = ( p - i Xi) ^ 


y 2 y4 

2 + V( + Tj. ) 

2(l-X-iZ) 4(l-X-iZ)2 


(4.4) 
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whcri^ 


a = real part of the phase refractive index 
Xj_- imaginary part of the phase refractive index 

Z = ^>/ m 


X = / co2 = 

IT ^ 


f 


IT 


ITfi- 


m SqW 


Y = 0 ) /(jj = eB/m w 
il 


Yj^ = Y cos 0= e Bg^/m u 
Y^ = Y sin 0= 


V = collision frequency of electrons 


fjj = plasma frequency 

f = operating frequency 
m = angular frequency 

e = charge of an electron (1.59x10“^® coulombs) 

B = magnetic induction vector 
m = mass of an electron (9.11x19"'-'^^ KgV 
IT = number density of free electrons 

= electric permittivity of free space (8. 854x10"^^ F/m) 

0 = angle between wave normal and geomagnetic 

field . 

In the absence of electron collision 


2X(1-X) 


T)' 


2-|liX)-Y^ ± ^(Y^, + 4Y^(1-X)^) 


(4.5) 


1 
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In the absence of magnetic field 

■A 


Tl 


2 


1 - 


X 

•i^iZ 


(4.6) 


leglecting both the magnetic field and electron 
collisions ; 

8111(h) 

= l-x = 1 (4.7) 

f" 

where N is electron density as a function of 
height . 

9 10 

4.1.2 Ionospheric Electron Density Model * t 

An analytical model, which represents the 
general features of the equatorial ionosphere 
(afternoon, equinox, sunspot maximum) was used to 
represent electron density contours. The model is 
an alpha Chapman layer with parameters which vary 
with geomagnetic latitude. It is represented by 

f2 ^ ^2 ^i(l-Z-g^) (4.8) 

K ” c 

where Z =(h - / H 

fj^ = the plasma frequency in IIHz 

f = the critical frequency in MHz 
c 

h = the height of ma.ximum electron density 

max 

in Em 

H = the scale height in Km. 



33 


f , h and H vary with the geoiragnetic latitude 

c ihclx 

in the f oil owing ways ; 

= 15 Ifflz for h< 100 Km 

fc = V'(50(X/8)2 exp(2-l X/8| ) + 40) for h>100Km 

h. = 350 Em for h< 100 km or 

nisix 

X >24° 

h = 430 + 80 cos( ~7 X) for h>100 Km and 

mo.x 4 

X < 24° 

H is d6;termined hy the constraint that f^^ = 2 MIfe 
at h = 100 Km 

X is the geomagnetic latitude in degrees. 

The profiles of ionospheric electron density model 
(BULGE) for four different latitudes are illustrated 
in Eigure 4«1* 

9 10 

4.1.3 Electron Density Perturbation Model ’ : 

The following model was used to represent 
the ionospheric electron density perturbation. It 
is an east-west toroidal irregularity with an 
elliptical cross section 

IT = 1Tq[ 1 +A (r,^^)] (A. 9) 


where 
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N = electron density of the ionosphere including 
the‘:ff.fects of: the torus 

Nq * electron density of the ionosphere neglecting 
toruseffect (specified hy any electron densjty 
model) 

^ (R.+H )(0-7r/24-X)cosp+(R-R -H )sinp * 

A (r, e) =: Co exp[-( -2—5 £—2 } 

A 

(R-R~HQ)cosp-(RQ+HQ)(9-Tr/2+X )sinp^2^ 

-( g ) j 

( 4 . 10 ) 


where 

= amplitude of the perturbation at the center 

R^ = the radius of the earth 
R, position in geomagnetic spherical polar 
coordinates 

A = the semiiiiajor axis of the ellipse in Km 

B 3S the semiminor axis of the ellipse in Em 

P - the tilt of an ellipse 
Hq * the height of the torus above ground 
X = geomagnetic latitude of torus. 

Throughout the remainder of this thesis ; the electron 
density perturbation will be referred to as the Torus. 

2 

4.1.4 Model of the Earth's Magnetic Field ; 

The earth's ma,gnetic field was represented by 
an earth centered dipole model. This is determined 
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from the following mathematical relationships; 

R 

I H I = Hq(— ( 1 + 3 cos" e r (4-11) 

6 = tan ^(2 cot 9) (4.12) 

where 

I H 1 = magnitude of earth's magnetic field 
Hq = surfaxe ma.gnitude of eaxth's mo.gnetic field 
at equator 

h “height above ground 
= Radius of earth 

0 = geomagnetic colatitude 

6 = magnetic dip angle 

4 . 2 Doppler Frequency Shift 

A radio wave propag?^ting through a medium 
whose index of refraction is changing with time 
experiences a. Doppler frequency shift. The 
instantaneous frequency shift A f» in cycles per 
second, is given by the time derivative of the 
phase shift ; 

Af = ~|^( A 9/2m) (4.13) 

The phase shift is obtained from 
A<P = 211 I AP 


(4.14) 
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where A P is ihe change in the phase 
(4*13) and ( 4 . 14) together give 

A f = ~ - A P) 
c dT^ ^ 


path. Equations 


(4.15) 


the negative sign indicates that a decrease in phase* 
path produces a positive frequency shift and vice versa. 
Prom Davies (1968) 


d( AP) 

di 


dl 


d 

dT 


/ h ds=^ ^ ^.ds (4.16) 
3T^ 


s s 

where s is the ray path, and is the real part of 
the phase refractive index. 


Substituting (4.16) into (4.15), we have 


“ c / ds (4.17) 

Neglecting the earth’s magnetic field and electron 
collisions, we have 

Ltt ~ (q _ 

3 T 3 

= • ®-| ( 4 . 18 ) 

Where K is a constant equal to 80.5x10^. The 
expression for A f may thus be written as 


Af = 


E , 9_N 

2^if 4 ^ ^ 


ds 


(4.19) 


As discussed in Section 2.3.4) the torus is thought 
to move in a North-South direction with an average 
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velocity of 160 m/sec. feking this into account and using 

dl'T' 

equation (4.9)» ~ may he conveniently derived from 


dT dx dx dT 


( 4 . 20 ) 


lA/here x is the longitudinal displacement of the torus 
relative to the geomagnetic pole, and X is the geomagnetic 
latitude of the center of the torus. 


12 

dx 


Since ~ equals - 


1 

r 


dN 


( 1 dg dx 
r '' dX dT 


hut H = n^[l +^ (r,9,(p)] 


( 4 . 21 ) 

(4.22) 


Therefore 

^ wr N d , ^ 

lii 2 — (^(r,e,<p)) ~ (4.23) 

dT r dX dT 

, dA (r,9,m) 

where ~ ^ is obtained hy differentiating 

equation (4.10) 


|_Mr,e,9) = . 2 i(^,e,,p) (m£0|-6 - I-Sfa-i) 

Cl X ° 


(4.24) 


where 


P = 


(Rq+Hq) (e-V2+^)cos P + (R-RQ-HQ)sin p 


y = 


(R-Rq-H^) cos P - (9-11/2 + X ) sin P 

— 


Thus equation (4.19) may he expressed as 


= 2^ / ft ■r'- 

s 


F P cos P Y sin p 


B^ 
(4.25) 


•)ds 
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since ds 


dt 


t is group path (independent variable) 


a(Af) 

dt 






( VHo) If A (r,e ,,p) 


B' 


or 


= (R +H ) (g-Cos. P Y sin 3 s ^ ^ « k 

dt ,2 ^ ^ A(r,G,9) — 

^ uu 


N dx 
o 


pp'f ^ 

(4.26) 

The above equation is in differential foric relative to the 


group path length t. To find Doppler frequency shift 
induced by the moving torus upon the radio wave, equation 
(4»26) is integrated simultaneously along the ray path 
with the Haselgrove equa.tions. 


4 • 3 Radio ¥ave Absorption 

The attenuation of radio waves in the ionosphere 
results from collision between electrons and other particles. 
The incremental absorption per tinit length (K) is expressed 
as 

jg; _ ^ (8.68x10^) dB/Km (4.27) 

where Xi is the imaginary part of the complex phase refrac- 
tive index, m is the angular frequency and c is the velocity 
of light. Neglecting the effect of the earth's magnetic 
field, Xi Diay be written as 



1 X Z ^ e^ 

1 + 2e^ m 0 ) p 

5 . 3x10'"^ N v/o) 


N V 


m" + 


(4.28) 


where the terms appearing in the equation (4.28) are 
as defined in Section 4.1.1. Therefore? 
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K = Q«046 U V 
^ 

The total absorption. A, 
evaluated from? 


(4.29) 

over the ray path, S, is thus 



0.046 




0)^ + 


V 


ds 


(4.30) 


Expressed in differential form relative to the group path length 
(t),the above equation, which is integrated simultaneously 
with the ray tracing equations, becomes? 


dA _ 
dt 


0.046 NV 

|ip,* OD ^ 


(4.31) 


4-4 Program Structure 

A three dimensional ray tracing program written 
by Jones (1966, 1968) in PORTRAIT language for the GDC 36 OO 
computer and modified by Sarabjit (1972) in FORTRAR' IV for 
an IBM 7044 computer has been used. It ha,s a main program 
followed by essential subroutines. When using this program 
the ionospheric mor' els, which define electron density, and 
the earth's magnetic field as a function of position in 
spa,ce, have been specified. The effect of electron 
collisions have been neglected. Each of these characteris~ 
tics of the ionosphere is defined by a. separate routine. 
There are several versions of the subroutine RIHDEX- These 
correspond to inclusion or exclusion of the effects of 
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earth, s magnetic field and collisions. Data required for 
initialization have been read into an array designated as 
the ¥ array. 

4 . 4.1 Accuracy of the Program; 

The numerical integra.tion subroutine has a built-in 

mechanism to check errors and adjust the integration step 

size a,ccordingly . The user specifies the accuracy he wants 

with W(42). To get a very a,ccurate but expensive ray trace 

—7 —8 

a small ¥ 42 (about 10 or 10 ) may be used, while 

for a cheap but approximate trace a large value of ¥ 42 
(10 ^ or even 10 ) is used, ¥ 42 is the maximum allowable 
relative error in any sin^e step for any of the equations 
being integrated. For the purpose of this work an accuracy 
of 10 ^ has been used. 

4.4.2 Coordinate Systems: 

The program uses two different spherical pola.r 
coordinate systems, namely, a geographic and a computa- 
tional coordinate system. Input data for the coordinates 
of the transmitter (¥ 14 and ¥ 16), and input data for the 
coordinates of the north pole of the computational coordi- 
nate system (¥ I 3 and ¥ 15) are entered in geographic 
coordinates. If ¥ 13 is equal to 0, ¥ 15 equal to 90*^, 
the two north poles superimpose and the two coordinate 
systems become the same. ¥hen the two coordinate systems 
do not coincide, the three type of ionospheric models 
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calculate electron density, the earth’s magnetic field, 
and collision frequency in terms of the computational 
coordinate system. In this work the dipole model of the 
earth’s magnetic field uses the axis of the computational 
coordinate system as the axis for the dipole field. Vdaen 
dipole model is used, the computational coordinate system 
is a geomagnetic coordinate system, and therefore electron 
density has been defined in geomagnetic coordinates. 
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Chapter 5 
RESULTS 

The comp'uterized. three dimensional ray tracing 
program, described previously, has "been used to study the 
ray path trajectories of the radio wave propagating 
through the ionosphere in HP hand of frequencies, effect 
of torus on the ray path trajectories, phase path/group 
path variations with elevation angle, and phase path 
oblique ionograms. The details and ressults of these 
predictions are presented in the following sections* 

5.1 Ray Path Tra.iect cries 

The results have been computed for the radio waves 
in HP band of frequencies (5*30 MHz) for all elevation 
angles of transmission, in presence of magnetic field 
and neglecting the effect of electron collisions. These 
computed results are given in Table 5.1- (In this table 
True Path is the path traversed by a radio wavewdtti free 
spacej^ Bufch the transmitter as well as receiver are 
located on ground, transmitter being located at Kanpur 
(26.5°H, 80° B) . The azimuth angle of transmission has 
been taken as 180° corresponding to the Horth- South 
propagation along constant geomagnetic longitude. The 
ray paths have been plotted using the IBM 1800 computer 
and associated plotter only for frequency range of 
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5-20 MHz at step of 5 Mz; and for elevation angle 30*^-90° 
at step of 10°. These plots are illustrated in Figures 
5.1 through 5.4. The ray paths plotted for the undistur- 
bed case exhibit typical features of ionospheric radio 
wave propagation, in that : 

- the ray paths are symmetrical about the point 
of reflection 

- the height of reflection increases with the increas- 
ing elevs-tion angle» consistent with- the electron 
density profile used 

- rays propagating vertically upwo.rds are returned 
directly to the transmitter 

- the rays bend towards the region of decreasing 
electron density; that is; rays bend towards the 
earth in the bottom ha.lf of the layer and away from 
the earth in the top half of the layer, etc. 

- a 20 MHz signal studied for 30-60 degrees at the 
step of 10 degree elevation angle penetrates through 
the ionosphere at 60 degree elevation angle. It 

is observed that the ground range decreases as the 
elevation angle is increased from 30 degree to 40 
degree. The rays which experience this effect are 
low angle rays. A 50 degree elevation angle ray 
is reflected back towards the ground at distance 
longer than the skip distance. This is higher . 





Fig. 5.2 Fre 





angle ray and is known as a pederson ray. For higher 
nngle rays, the ground range increases with increase 
in the elevation angle. The mechanism of pederson ray 
runs as given helow. 

The refractive index decreases as the wave 
penetrates "into the regions of greater electron 
density (equation (4*7)) and the angle of refraction 
increases correspondingly (equa,tion ( 2.6)). It is 
seen from Bulge profile of Figure 4*1 that dU/dh, the 
rate of change of electron density with height, decreases 
until it hecomes zero at the height of maximum electron 
density and beyond this height it becomes negative 
corresponding to the penetration condition. The rays 
reflected back towards the ground ffom the height of 
maximum ionization experience small bending as dU/dh 
is very small and return back to the ground at distance 
larger than the skip distance. These rays are high 
angle rays and are known as pederson rays. 

— the critical frequency at Eanpur is 15*4 MHz. 

The radio waves transmitted at frequency higher 
than 15.4 MHz penetrate through the ionosphere at some 
elevation angle less than 90 . In these ca.ses there 
will be both low angle as well as high angle rays. 
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5 • 2 Effect of Perturbation on HF and VHF Propagation 
To study the effect of torus on HE radio wave 
a 20 MHz signal has "been considered for elevation angle 
30—60 Deg at step of 10 Deg. The torus is located in 
space such that radio wave pa,sses through the center of 
the torus. The torus to background ionization density 
has been taken as 1.1. The semimajor and semiminor 
axes of an ellipse have been taken as 150 Km and 50 Km 
based on the experimental observations. The torus is 
tilted down to the South at an angle of 15 degree. 

Both the transmitter and receiver are located on the 
ground, transmitter being located at Kanpur (26.5 H, 
80°E). To observe the effect on VHP signal, a 41 
MHz signal transmission, from satellite located at the 
height of 1000 Km above ground, has been considered. 

The receiver in this case is located on ground. The 
computed results are shown in Table 5.2 and the plots 
for HF signal are illustrated. in Figures 5.5 and 5.6. 

The effect of the presence of the torus on the radio 
wave path are listed below. 

- Ray path trajectories are significantly affected 
in the vicinity of the torus, the rays whose wave 
reversal occur near the center of the torus are the 
most significantly affected. 
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The ray path symmetry ahout the point of reflection 
is not significantly disturbed for the intensity of 
perturbation used in this work. 

There is significant change in various path lengths 
and skip distance. Por example, for a 20 JVIHz ray 
with an elevation angle of 50°, the phase path and ground 
range can change by as much as 100 Km depending upon the 
location of the torus. 

As the torus is moved away from the transmitter towards 
the point of reflection, where the ray path between 
ground and the ionosphere experiences the presence 
of the torus, the various path lengths, ground range 
and maximum height of reflection of radio wave decreases 
and this decrease in various parameters is maximum when 
the torus is located near the point of reflection. If 
the transmitter torus separation is still increased 
the radio wave now experiences the presence of the 
torus along its return path to the ground. The 
effect of the torus now decreases and finally the 
radio wave path in the presence of the torus is 
almost same as an undisturbed path. 

The effect of torus on VHP radio signal transmitted from 
■^satellite at 41 MHz is insignificant as can be seen 
from Table 5.5. (In this table straight line path is the 
path joining the transmitter and receiver by a straight 
line and Pifferential Phase Path is the difference between 
straight line path and the Phase Path). However, this 
small change in Phase Path and Straight line Path is 
significant in radio astronomy. 
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5 • 5 Variation of Phase Path/G-roup Path •with Elevation 

6 

An^le 

In this case the effect of earth’s magnetic field 
as well as collision have been neglected a,nd computed 
results for phase path, group path, and ground range for 
various eleva,tion angle and for frequencies below and 
above critical frequency are shown in Table 5«4. The 
curves, depicting the effect of elevation angle on phase 
path/group path and ground ro.nge for a fized frequencj’’ 
below and above critical frequency, have been shown in 
Figure 5.7 and Figure 5.8. It is observed that 

- For ’f (13 MHz) less than vertical critical frequency 

f (15.4 MHz) phase path (P) tends to twice the phase 
c 

height as elevation angle approaches 90 and D approaches 
zero. 

- the tangent to the curve (see Figure 5.7) is horizontal 
at D equal to zero. For frequency ‘f’ (18 MHz) there 
is a minimum value of ground range (the skip distance ) 
at which phase path is. also minimum. 

- The low angle and high angle branches of the curve 
meet in a cusp as illustrated in Figure 5.7. 

~ A plot of effective (group) p.ath P' ag9.inst ground 
distance D shows a ’nose' near the skip distance 
(see Figure 5.8). The behaviour of the (D ,P) curve 
is thus quite different from that of the (B,P' ) curve. 


2000 
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Fig. 5. 7 Variation of D and F’ with elevation angle. Bulge, f=13MHz 
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- The slope of (D,P) curve does not undergo any large 
changes near the cusp, where as the slope of (D,P' ) 
curve changes rapidly round the nose. 

- These observations are in agreement to those found 
by Davies (1968) and Gething (1974) for a single 
parabolic layer. 

5*4 Ph ase Pa,th Obliqu e lonogram^ 

Por a. given ground distance, phase path/group path 
oblique ionogram is a plot of phase path/group path 
against , frequency for different elevation angle. Por 
frequency of operation below critical frequency there 
is one set of elevation angle and frequency for which 
a radio wave is received between two fixed terminal 
stations, whereas for frequency of operation greater than 
critical frequency there will two sets of elevation angle 
and frequency of transmission, for the same communication 
links. The ray tracing program developed by Jones (1966, 
1968) does not incorporate automatic homing technique 
to select a set of elevation angle and frequency of trans- 
mission for a given ground range. However, this may be 
done by step-by-step method by getting results computed 
for the nearest Kilometer accuracy in ground range. In 
this way the user has to pay an exorbitant cost of 
computation. In view of all these difficulties phase 
path oblique ionogram could not be completed. The results 
obtained for some sets of elevation angle and frequency 
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of transmission for two different groiind ranges 500 Km and 
250 Km, are given in Table 5 . 5 . While computing these 
results efforts have been mode to select a set of elevation 
angle and frequency to get computed ground range nearest 
to the specified value (250 Km and 500 Em) for a given 


communicant ion link. 
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TABLE 5.1 COICPUTATIQKF OP GROUFD RANGE, PHASE PATH AKI) GROUP 
PATH FOR VARIOUS FREQUENCIES AND ELEVATION ANGLES 


apogee 

Deg. EM 

GROUND 

distance 

KM 

GROUP 

PATH 

KM 

PHASE 

PATH 

TRUE 

PATH 

KM 


FREQUENCY = 5 

MHz 




0 

65.8538 

2107.3824 

2125.97180 

2117.05185 

2121.11102 

10 

78.1543 

935.8076 

965.67976 

949.42618 

956.34807 

20 

94.7220 

607. 1057 

659.43687 

651.22401 

643.82150 

30 

109.1143 

465.3888 

551.60091 

502.56491 

523.69719 

40 

120.9260 

561.9711 

486.55055 

417.27026 

445.95750 

50 

129.9931 

277.9950 

444.82224 

357.26753 

391.49256 

60 

156.7102 

204.1852 

417.86274 

514.15054 

351. 31815 

70 

141.2751 

133.5909 

599.74964 

285.52249 

320,01525 

80 

143.7159 

66.1597 

588.96741 

264.67869 

298.09356 

90 

144.2511 

5.2528 

384.49092 

258,00597 

288.43306 


FREQUENCY = 10 MHz 




0 

91.9760 

2477.9446 

2508. 25151 

2495.08606 

2501.36630 

10 

110.2057 

1513.7575 

1560.96170 

1334.53815 

1347.08311 

20 

135.7507 

908.2615 

989.20760 

941.08037 

963,56448 

30 

156.5410 

689.1535 

814. 75109 

739.93250 

773.91609 

40 

173.8873 

540.5857 

720.75994 

615.24419 

661*21008 

50 

188.2905 

422.1084 

666. 24961 

528.08496 

584.81467 

60 

199.7578 

516.8574 

653.82719 

465.86781 

528.13512 

70 

207. 9897 

215.1367 • 

615.84692 

416.57449 

483.02482 

80 

212.5782 

111.8975 

600.87536 

585.05031 

446.85455 

90 

215.3033 

12.8560 

593.30497 

571.56338 

426.90771 


(Continued) 
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PREQUEFCY 15 MHz 


0 

140.8047 

3414.3648 

3499.77554 

3444.06110 

3470.98087 

10 

143.0990 

1814.7361 

1883.41554 

1843.13948 

1862 . 52307 

20 

171.7526 

1199.9516 

1307.17952 

1241.26587 

1272.48247 

30 

199.1646 

922.9800 

1088.70802 

983.37553 

1031.87770 

40 

224.9908 

748.4447 

991.99008 

833.55411 

903.44699 

50 

249.7309 

613.7998 

957.91555 

730.99762 

825.06342 

60 

272.9041 

492.5008 

965.83663 

652.96514 

771.56580 

70 

293*8047 

367.5669 

1006.39594 

590.03139 

728.49223 

80 

509.5821 

217.9155 

1064.44603 

538.35811 

681.05576 

90 

314.6039 

39.4984 

1095.98782 

507.40891 

632.39747 

EREQUENCY 20 MHz 




0 

23^.8679 

3310.0286 

3486.55981 

3363.74738 

3422-. 7111 5 

10 

177.0793 

2829.9429 

2927.12329 

2869.42007 

2897.38443 

MIN. 

.DIS 23.8780 RAY MADE A CLOSEST 

APPROACH TO 

REC HT. 

20 

209.0555 

1590.0499 

1728.52063 

1639.05898 

1681.71149 

50 

247.6866 

1261.6192 

1480.40762 

1326,90823 

1389.08121 

40 

293.8677 

1139 . 0911 ® 1489.24504 

1217.35187 

1337.93274 

50 

377.4663 

1477.8815 

2153.32190 

1498.10825 

1770,59962 

60 

413.8076 

320.6792 

676.25000 

445.61528 

537.99447 


RAY PENETRATES 


(Continued ) 
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FEEQUEKCY = 25 Mz 

0 541.1835 3666.8344 4034.52170 3724.58898 3869.07370 

10 220.1012 3882.7427 4063.91681 3951.06064 4006.01123 

20 255.3998 2789.0248 2974.68396 2846.97586 2908.17441 

30 323.6197 2578.1222 2900.70850 2635.19073 2759.16971 

40 455.8933 752.5902 1045.00000 814.95699 917.80948 

R/VY PENETRATES 
ERBQUBNCY = 30 MHz 

0 471.2875 2865.6852 3063.75000 2923.62338 2990.33286 


RAY PBNETRiiTES 



TABLE 5.2; EFFECT OF PERTURBATION ON G-ROUIID RANGE, 

PHASE Path anl group path 
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Frequency 20 MHz 


Eleva- 

tion 

Angle 

Deg. 

APOGEE 

KM 

GROUND 

DISTANCE 

KM 

GROUP PHASE 

PATH PATH 

KM KM 


TRUE 

PATH 

KM 

Mag. Field, Extra 

ordinary, No col.. No perturbation 



30 

247.6866 

1261.6192 

1480.40762 1326.96823 

1398.08121 

40 

293.9795 

1139.0910 

1489.24480 1217.35162 

1337.93250 

50 

377.4663 

1477.8815 

2153.32190 1498.10825 

1770.59962 

60 

413.8076 

320.6792 

676.25000 445.61528 

5 

37.99447 

RAY PENETRATES 





No Field, No col.. No perturbation 



50 

423.5107 

3038.6431 

3929.72089 2971.42166 

3587.24152 

Torus- 
s ep era 
tion 

KM 

Tx Height 

of torus APOGEE 
above 

ground ^M 

KM 

GROUND GROUP PHASE 

DISTANCE 1’'^™ PATH 

KM KM . KM 

TRUE 

PATH 

KM 

Mag. field. No. col., Perturbation 10 percent. Elevation 50 Deg. 

220 

250 

365.772 

1355.120 1985.595 1392. 

972 

1640.192 

440 

350 

357.718 

1294.831 1915.536 1338. 

352 

1576.868 

660 

380 

374.507 

1239.316 1890.452 1294. 

672 

1539.879 

880 

330 

377.487 

1441.000 2126.426 1460. 

960 

1736.718 

1100 

240 

377.466 

1477.324 2156.245 1495. 

650 

rr 

1770.486 


No field, No col. 7 Perturbation 10 percent, 50 Deg. 

455.943 1906.654 2684.738 1782.139 2166.860 


660 


380 



TiSLl 5.3: IONOSPHERIC EFFECTS ON SATEIilTE SIGNALS 
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TABLE 5.4: DATA FOR PHASE PATH VERSUS GROUND RANGE 
CHARACTERISTICS 


El eva- 
tion 
Angle 
Deg. 

PHASE 

PATH 

KM 

GROUP 

PATH 

KM 

GROUND 

RANGE 

KM 

APOGEE 

KM 

No field, no col. 
FREQUENCY = 13 MHz 




10 

1660,11203 

1694.82843 

1634.3637 

132.5854 

20 

1137.69321 

1195.68356 

1099.0154 

160.0689 

30 

897.15826 

988.60684 

839.3936 

185 .-0800 

40 

752.74t68 

886.65970 

669.7663 

207.5634 

50 

651.60728 

855.86947 

535.2146 

227.6503 

60 

574.59824 

814. 52509 

412.6176 

245.2708 

70 

514.53995 

809.48001 

289.7041 

259.1954 

80 

472.18711 

810.70087 

160.2658 

268.1114 

90 

451.51850 

807.73646 

24.7307 

270.3928 

No field, no col. 
FREQUENCY = 18 MHz 




10 

2446.0235 

2480.89893 

2530.82797 


20 

1440.4579 

1487.01300 1566.08089 


30 

1126.6125 

1191.18713 1522.78902 


40 

961.0701 

1045.89752 

1261.74823 


50 

884.1586 

984.52897 1352.54863 


60 

2556.8380 

2374.71051 3712.71152 
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TABLE 5.5- DATA" BOB PHASE PATH OBLIQUE IOHOGRA.M 


Elevation 
Angle Deg. 

GROUND ~ 
RANGE 

IM 

“Phase 

PATH 

EM 

GROUP 

PATH 

KM 

EREQUENCY 

MHz 


No 

field, no 

col . 



10 


501.5768 

504.63325 

512.6O515 

0.91 

20 


500.1307 

518.43672 

539.68216 

2.73 

30 


500.4471 

539.83060 

589.45547 

5.36 

40 


499.8812 

571.41452 

664.67851 

8.55 

50 


500.5769 

615.05274 

783.21387 

12.10 

60 


500.8456 

662.14441 

979.32755 

14.90 

70 


500.5505 

693.20570 

1342.77449 

16.06 

10 


250.3083 

249.67565 

255.63869 

0.285 

20 


249.4446 

254.05931 

267.43227 

0.485 

30 


250.1973 

267.93432 

292.05846 

0.990 

40 


250.6223 

291.14130 

331.54773 

2.070 

50 


249.7085 

323.50725 

395.69491 

3.800 

60 


249.8515 

380. 13170 

503.01260 

6.900 

70 


249.2546 

465.80571 

703.61874 : 

11.500 

80 


250.2635 

565.98330 

1209.69995 

15.230 
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Chapter 6 
COFCIUSIOKS 


A three dimensional ray tracing computer program 
has been used to study the following effects of ionos- 
pheric propagation: 

(a) Effect of elevation an^e and frequency on HE 
radio waves through the ionosphere. 

(b) Effects of ionospheric irregularities on ground- 
to-ground HE rays. 

(c) Effects of ionospheric irregularities on Satellite- 
to-ground VHE rays. 

(d) Phase path characteristics for oblique angle rays*. 

In addition computations were made to. study the Doppler 
shift produced by motion of ionospheric irregularities. 

The following conclusions may be drawn from the 
results of this simulation. 

1. The ray path trajectories through an undisturbed 
ionosphere are symmetrical with respect to the point 
of reflection. 

2. Ray path trajectories are significantly affected 
in the vicinity of the electron density irregularity. 
The irregularity model chosen was that of a torus; 
the rays whose wave reversals occur near the center of 
the torus are the most severely affected. 
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3. As the torus transiPitter separation increases, the 
effect Of torus on ground range and various path lengths 
increases and the effect is maximum when the torus is 
located near wo.ve reversal point. It is observed ths,t 
focussing -of high angle radio waves with low angle 

radio ways, in the presence of torus, may result in 
increase or decrease in the amplitude of received signal. 

4. The effect of the torus on VHP signal transmitted 
from satellite is insignificant. 

5 . The Phase Path versus G-round Eange chara,ctcristic 
for various elevation angles has a zero slope for elevation 
angle 90°. 

6. For frequency of operation greater than vertical 
critical frequency the low angle and high angle branches 
of the curve, phase path (P) versus ground range (D), 
meet in a cusp. A plot of group path (P' ) against ground 
range (D) shows a 'nose' near the skip distance. 

7. The slope of the (P,D) curve does not undergo any 
large changes near the cusp, whereas that of (P',D) curve 
changes rapidly round the nose. 

8. There was no noticeable Doppler shift for sr^tellite- 
to-ground waves at 41 MHz when the torus velocity was taken 
as 160 m/sec. 

Various other topics such as phase distortion of 
UHP signal, absorption, Doppler shift, focussing and 
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d ef ocussing of radio waves d.ue to ionospheric irregular- 
ities can he studied using ray tracing techniques. For 
complete ionospheric probing a reader should also consider 
the effect of lower atmosphere (troposphere) particularly 
when considering frequencies in the microwave band. The 
communication from satellite to another satellite can be 
studied in which case Doppler effect will be of 
importance . 
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